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Abstract 

Magnetostatic normal modes are investigated here in elongated rods. The dipolar field resulting 
from the dipole-dipole interactions is calculated numerically in points of the axis connecting oppo- 
site rod face centers {central axis) by collecting individual contributions to this field coming from 
each of the atomic planes perpendicular to the central axis. The applied magnetic field is assumed 
to be oriented along the central axis, and the magnetization to be uniform throughout the sample. 
The frequency spectrum of magnetostatic waves propagating in the direction of the applied field is 
found numerically by solving the Landau-Lifshith equation of motion with the spatially nonhomo- 
geneous dipolar field taken into account; the mode amplitude profiles are depicted as well. While 
energetically highest modes have bulk-extended character, the modes forming the lower part of the 
spectrum are localized in the subsurface region (bulk-dead modes). Between these two mode types, 
magnetostatic modes of a new type (comb modes) are found to occur, characterized by two clearly 
discernible regions: a zone of fast amplitude oscillations inside the rod, and narrow slow-oscillation 
regions at the borders. Absorbing virtually no energy from an applied alternating field, comb 
modes will have no significant contribution to the magnetic noise. 

PACS numbers: 75.30.Ds, 75.40.Gb, 75.75.+a 
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I. INTRODUCTION 



Modern magnetic storage media require the use of very small magnetic elements (see e.g. 
Q]). When the recording density approaches 100 Gb/in 2 the bit size is around 0.2//mx0.1//m, 
and similar size is required for the read heads. On the other hand, when the bit size 
becomes smaller, the sensitivity of the read heads should increase in order to detect the 
weaker magnetic field produced by the bit. In small nonellipsoidal magnetic particles the 
internal field nonhomogeneity, due to the particle shape and produced by long-range dipolar 
interactions, has to be taken into account. It is known that dipolar size and shape effect 
present in the magnetostatic wave spectra of such particles result in additional power losses^ 
very important in device applications. Therefore, one way to improve the quality of the read 
heads is to suppress their magnetic noise. In this paper we show that the magnetostatic 
mode spectrum of the elongated nanorods reveal specific range of frequencies when new type 
of excitations (we shall call them " comb modes" ) are accessible and when they are excited 
by the AC magnetic field the absorbed power is negligible. Therefore, in this particular 
range of frequencies, magnetostatic modes will not contribute to the level of the magnetic 
noise. 

The dynamical properties of small magnetic elements have been actively investigated re- 
cently in studies using different methods and approaches (see our recent papers and 
references therein). These different methods give good indications of the localized magneto- 
static mode existence. This paper is focused on localization properties of pure magnetostatic 
modes in elongated rods (having a nanometric square cross-section). It should be emphasized 
that the localization effects we are going to discuss here can be revealed only in magnetic 
systems in which dipolar interactions are dominant, and for this reason, exchange interac- 
tions are omitted here from the very beginning. If exchange interactions were taken into 
account in our investigations, the mode localization properties in question would disappear 
completely 

In our present theoretical approach the studied dipolar system is regarded as a set of 
discrete magnetic dipoles regularly arranged in a crystalline lattice. The dipolar energy 
is calculated by collecting contributions from each dipolar lattice plane parallel to the rod 
base (see Fig. [TJ. As the summing is performed over all the dipoles within a given dipolar 
plane, no approximation is involved in our dipolar energy evaluation. The scheme of our 
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calculations has already been presented in our former paper [4j , restricted to cubic systems. 

II. A FINITE SYSTEM OF MAGNETIC DIPOLES IN PLANAR ARRANGE- 
MENT 

We shall consider a dipolar system consisting of magnetic moments arranged regularly 
in sites r of a simple cubic crystal lattice and coupled by the dipole-dipole interactions. 
The system is assumed to form a rectangular prism with square base (see Figure^). Let 
the prism base determine the (x, y)-plane of a Cartesian reference system, with the z-axis 
perpendicular to this plane. The reference point (0, 0, 0) shall be placed in the central site 
of the prism bottom. 

Let us calculate magnetic field hs "produced" by all the prism dipoles in a site indi- 
cated by internal vector R. According to the classical formula (obtained using the linear 
approximation field can be expressed as follows (in the SI units): 

r 1 3(r-R)U,-(r-R))-fI,\r-R\ 2 

^ = ^ E \r—R> ' (1) 

r+R I' ""-I 

the above sum involving all the sites except the reference point (i.e. the site with position 
vector f=R). The lattice planes parallel to the prism base shall be numbered with index 
n G (0, iV — 1) (see Figured), and the sites within each plane indexed with vector fn = 
a\pi + qj], defined by integers p, q G (—L,L). This means that the position of the sites in 
which the magnetic moments are located, indicated by vector r, shall be defined by a set of 
three integers, (p,q,n): 

r = [rj|, an] = a\p, q, n], p,qe(—L,+L) and n G (0, iV — 1), (2) 

a denoting the lattice constant. Thus, the considered prism contains N(2L + l) 2 magnetic 
moments. Below we shall focus on the magnetic field on the z-axis only, assuming its 
direction to be solely allowed for magnetic wave propagation. Hence, we put R = a[0, 0, n'), 
where n' G (0, N — 1), and re-index the dipole field: h n > = h^. Additionally, we shall assume 
that all the magnetic moments within a single planen are identical, i.e.: fl n = fl[ Pt q,n], for any 
p and q. It is convenient to introduce here the notion of magnetization, a phenomenological 
quantity, which in the considered case of simple cubic lattice can be defined as follows: 
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FIG. 1: (a) The rod sample considered here; the applied field, Hq, sets the magnetic moments in 
the direction perpendicular to the rod base (i.e. along the z-axis). The rod 'thickness' is (N — l)a, 
and the square base side width is 2La (a denoting the lattice constant). The magnetostatic waves 
are assumed to propagate along the z-axis (i.e. in the direction of the applied field), (b) The 
planar rod model used in our calculations. 



M n = fl n / a 3 . Then, 



becomes: 
1 



Air 



(3) 



where we have introduced a matrix D whose elements, D n y, are defined as follows: 

,\(V 2 + q 2 ) - (n- n'f 



(4) 



p,q \<p2 _|_ q 2 _|_ ( n _ n r 

It should be remembered that site (0,0, n') is excluded from the sums appearing in the 
equation (jlj). Equation (J3j), defining elements D n<n < y can be interpreted as the definition of 
a dipolar matrix D composed of these elements; as we shall see this matrix shall play an 
important role in deducing properties of magnetic modes. As a consequence of (@J) matrix 
D is symmetric, i.e. D n ^ n i = D n ^ n . Let us introduce a new variable, defined as follows: 
5 = n — n'\ 5 measures the distance between planes n and n' . Consequently, we can write: 

>,W + q 2 )-5 2 



D 



n,n' 



Ds = J2' 



P,g [p 2 + q 2 + S 2 } 2 

Note that the considered system consists of a finite number of planes (index n taking 



(5) 



values n 



0,1,2, 



N — 1), so the set of values available to 5 depends on the reference 



plane, n', with respect to which the distance is measured. However, the following condition 
must always be satisfied: 

0<ri + 6<N -1. (6) 

A sum over all the system planes shall be in use below; this sum shall be denoted as: J2s n ' ■> 
superscript n' indicating that the summing is performed on the planes neighbouring with 
n' , including plane n' . Therefore, 5 takes the following values: 5 = 0, ±1, ±2 ... its lower and 
upper limits being determined by condition (0). 

Up to now, the direction of dipole arrangement has not had much importance in our 
reasoning. Now we shall consider the case with dipoles arranged along the z-axis only. 



III. DYNAMICS OF THE DIPOLAR SYSTEM 



In this paragraph we shall consider a magnetic prism placed in a static magnetic field, 
H , applied along the z-axis (Figure HJ). Field H is assumed to be strong enough to arrange 
all the magnetic moments along the z-axis. Then, the magnetization vector can be regarded 
as a superposition of two components: static (parallel to the z-axis) and dynamic (lying in 
the (x, y)-plane): 

M n = M s k + m n ; (7) 

Ms is the static magnetization, assumed to be homogeneous throughout the sample, and 
vector m denotes the dynamic magnetization, perpendicular to Ms- Similarly, the dipole 
field, h n i, can be resolved into two components: static, h s n , (parallel to the z-axis), and 
dynamic, h^, (lying in the (x, y)-plane): 

h n , = h s n , + h d n ,. (8) 

These two components of the dipole field can be easily found from (JHJ). By replacing the 
third component of the magnetization vector with the static magnetization (i.e. by putting 
M* = Ms), and the two other components, M* and M%, with the respective components of 
the dynamic magnetization, m x n and m y n , the following formulae are obtained: 



- E Dn,n> 



M s k, h d n , = ^-Y, D n,n>m w (9) 
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element D n ^ n i being denned by (HJ). 

The magnetic moment dynamics is described by the phenomenological Landau-Lifshitz 
equation (LL): 



dM 



R 



dt 



1^ M R >< H eff,R 



(10) 



eff r denoting the effective magnetic field acting on the magnetic moment in site R. This 
effective field is a superposition of two terms only: the applied field, Hq, and the field h n i 
produced by the magnetic dipole system: 



H e ff,R = H e ff, n ' = H + h n > = (H + h s n ,)k + h d n ,\ 

the above-introduced dipole field components (static and dynamic) being defined by 
The LL equation becomes: 
drh n i 



(11) 



dt 



-THo (M s k + m n ) x {{H Q + h s n ,) k + hi) 



(12) 



We shall solve it using the linear approximation, i.e. neglecting all the terms with m 



squared. Assuming the standard harmonic time-dependence of the solutions: m, 
(fT2"j) becomes: 



/ ~ e 



— iLom n > = 



M s k x h d n , + 7?v x (H + h s n ,)k 



or, using the properties of vector product: 

- iujm n i = 7/i A; x M s h% - m n i (H + h s n ,) 



(13) 



(14) 



Through replacing the dipole field dynamic and static components with their explicit ex- 
pressions © we obtain: 



iuim n i = 7/ioA; x 



X 

M s-r- D n,n>m n - m n , 



(ft-^EA^)], (15) 



or, after bilateral multiplication by — 47r(7/i Ms) 1 : 



iftm„,i = k x 



>m r , 



(16) 



Q and Qh denoting the reduced frequency and the reduced field, respectively, defined as 
follows: 



7/i M s 



and Qh 



AtiHq 

~m7 



(17) 



Two complex variables are now introduced for convenience: = ± im y n \ with these 
new variables, (fTHj) splits into two independent identical scalar equations for m+ and m~; 
this means we are dealing with magnetostatic waves polarized circularly. Therefore, it is 
enough to consider only one of these two equations, e.g. that for m+: 

Qm+ = m+ (n H ~ D n ,A - ]T D n , n ,m+ (18) 

\ n / n 

The above equation can be rewritten as follows: 

Qm+ = m+H(n') - £ D s m+ ±s , (19) 

5=1,2... 

where we introduced the following abbreviation denoting the local field: 

if(n') = n H -D -2Y, n 'D 5 = tt H + ^ d + W n ,, (20) 

5=0 

with n d = —D and Q^' = ~ 2 J2s=o n ' D§ meaning the {reduced) contributions to the local 
field coming, respectively, from the dynamical and statical parts of the demagnetizing field. 
Note that, the eigenvalues Q (being reduced frequencies) correspond to magnetostatic waves 
propagating in the direction of the applied field, i.e. along the central axis shown in Figure 
□ 

In the remaining part of our work we will be considering only the rod-shaped samples, i.e. 
starting from this point we always assume 2L < N — 1. We also assume particular values for 
fioHo = 0.2T and Ms = 0.139 • 10 6 AmT x (YIG magnetization) with resulting value for the 
reduced field Qh = 14.374. However, we have to emphasis that selection of this particular 
value for Qh is not essential for results to be presented in subsequent sections of this work, 
since the distribution of eigenvalues Q and profiles of modes associated with them are not 
sensitive to the choice of particular Qh value: the particular value of Qh only sets the whole 
spectrum in a given frequency region and if Qh changes the whole spectrum is shifted to 
another region, but the relative distribution of mode eigenfrequencies remains unchanged. 

IV. MODE FREQUENCIES AND AMPLITUDE PROFILES 

We shall investigate magnetostatic excitations in a rod of the plane size 2La = 20a (a 
being the lattice constant). The rod consists of ./V planes normal to the z-axis and numbered 
with index n, ranging from n=0 to n = N — 1, as indicated in Figure H The effective dipole 
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FIG. 2: The discretized magnetostatic mode frequencies vs. the rod length, (N — l)a; m indicates 
the mode number. 

field calculation procedure applied in the previous paragraphs allows to find the field in 
z-axis points only, i.e. along the rod central axis passing through opposite rod face centers; 
this is the idea of the approximation used throughout this study, and henceforth referred to 
as central-axis approximation. 

The whole sample is assumed to be magnetized uniformly (the corresponding magneti- 
zation value being Ms) and in the direction of the external field, applied along the z-axis. 
The z-axis direction shall be also the only one allowed for propagation of the magnetostatic 
waves studied in this paper. With these assumptions, the problem of motion - to be solved 
on the basis of equation fllOl) - reduces to a single dimension in the space of variable n; the 
domain of the investigated motion is the interval n G (0, N — 1), between two opposite rod 
face centers. 

Figure |2 presents the discrete spectra of numerically calculated magnetostatic mode fre- 
quencies in a rod of variable length; the spectrum evolution with increasing N is visualized 
by the depicted frequency branches, each corresponding to one mode of a fixed number m. 
The plot shows clearly that the frequencies stabilize after reaching some critical length; we 
note that lower modes stabilize at lower "critical" values of length N than higher modes. 

In order to discuss in detail the properties of magnetostatic mode spectrum in rods, let's 
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consider a sample with dimensions 20a x 20a x 200a, i. e. a square rod with length exceeding 
ten times the base side. The plot shown in Fig. El indicates that the mode spectrum 
corresponding to this aspect ratio is stable enough to assume that no significant changes 
would result from further rod elongation. Fig. shows eigenmode profiles numerically 
calculated in the considered rod. Though showing only selected modes, the depicted mode 
profile sequence provides full information on amplitude variation with mode frequency. The 
features of the mode profiles shown in Fig. allow to classify all the modes obtained into 
four qualitatively distinct groups: (a) slow-oscillating bulk modes, comprising ca. 20 modes 
with highest frequency values; the number of oscillations grows as mode frequency decreases; 
(b) modes m = 170 to m — 90, with different character of amplitude profile oscillations in 
two distinct regions: a central region with fast-oscillating amplitude, and outer regions with 
slow-oscillating amplitude; (c) modes m = 3 to m = 80, whose amplitudes are completely 
suppressed in a central region, communicating on each side with a fast-oscillation zone 
followed by a slow-oscillation region; (d) two surface-localized modes, m — 1 and m = 2, 
with amplitudes concentrated almost exclusively on the sample surface. For the purposes 
of further analysis, the mode amplitude absolute value profiles are also depicted in Fig. Eb- 
The mode profiles shown in Fig. ®L,b allow to identify the modes in groups (a) and (d) as 
backward bulk extended modes and surface localized modes, respectively; both these mode 
types are well-known from literature. The remaining two groups, to our best knowledge 
not yet reported in literature, represent a new quality, and therefore need to be discussed 
in detail. Modes from group (b) - henceforth referred to as comb modes, because of the 
characteristic 'comb' profile of their amplitudes in a central region - are characterized by 
two wave numbers: high in the central region, low in the outer regions. A characteristic 
feature of modes from group (c) is complete amplitude vanishing in a central 'zero region'; 
since this zero region broadens as mode frequency decreases, the non-zero amplitudes in 
the adjacent outer regions are thinning down as they descend. For this reason, modes from 
group (c) will be henceforth referred to as bulk-dead modes. As we shall see further into this 
paper, comb and bulk-dead modes appear as a result of rod elongation along a distinguishing 
axis. 
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FIG. 3: Numerically calculated magnetostatic mode profiles in a rod of length 200a depicted along 
the central axis (indicated in Figure ^) showing separately: (a) the dynamical magnetization m + 
relative values, and (b) the corresponding absolute values, \m + \. 

V. INTERNAL DIPOLAR FIELD NONHOMOGENEITY VS. COMB AND 
BULK-DEAD MODES 



Now let's superimpose the calculated magnetostatic mode frequency spectrum and the 
profile of local field H(n) plotted along the rod axis (Fig. EJ). A strong spatial inhomogeneity 
of the local field is apparent from the calculated field values; two regions, qualitatively 
distinct in terms of inhomogeneity, can be distinguished: outer regions, with steep profile 
bias indicating high local field gradient, and a central region, in which the profile bias is 
relatively mild. Note also that the average value of local field, H av , can be regarded as the 
border line between these two regions. The mode profiles, ordered by growing frequency, are 
superimposed on the local field curve, with energy scale (on the vertical axis) maintained 
(i.e. common for the mode frequencies H m and the local field profile H(n)). An interesting 
correlation occurs between the mode type and the position of the mode frequency with 
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respect to the local field profile. All the modes with frequency values above the local field 
maximum (H max ) are bulk modes. Those with frequency values between H max and H av 
are comb modes, and those with frequency values below H av are bulk-dead modes. Note 
that the fast-oscillation segment in comb mode profiles coincides (approximately) with the 
low-gradient segment in the local field profile; at the same time, the borders of the zero 
region in bulk-dead mode profiles follow the local field profile in its high-gradient (outer) 
parts. Also, note a mode density peak in the immediate vicinity of the average local field 
value H av ; as much as 150 mode frequencies, representing approximately 3/4 of their total 
number, are concentrated around this point! Further interesting conclusions can be drawn 
from superimposition of the complete mode frequency spectrum on the local field profile 
(Fig. Ek). The first thing to be observed is a double degeneration of bulk-dead modes (i.e. 
the modes below H av ); the meaning of this degeneration of mode energy levels is explained 
by Fig. Eb, showing profiles of all bulk-dead modes with their spatial symmetry specified. 
Clearly, the degeneration of bulk-dead mode energy levels is associated with that of the 
corresponding mode amplitudes: in each degenerated pair, one state is symmetric and the 
other antisymmetric, both showing (see Fig. Efc>) identical amplitude absolute value profiles. 
Mode degeneration disappears in the comb mode frequency range (which is understandable, 
as symmetric states in this group have non-zero amplitudes in the rod center, and thus 
amplitude degeneration is impossible). Returning to Fig. an interesting thing to note 
is the linear character of comb mode frequency variation with mode number m; the range 
of comb mode dispersion linearity falls in the immediate vicinity of H av , the value at which 
- as we have already anticipated on the basis of Fig. |U - mode density peaks (see the 
supplementary plot attached to the main one on its right side). 



VI. MAGNETOSTATIC MODE RESONANCE SPECTRUM AND OUTLOOKS 

The presence of comb and bulk-dead modes in a rod spectrum results in a very interesting 
resonance effect, illustrated in Fig. The top part shows the SWR spectrum calculated for 
a 20a x 20a x 200a rod; the modes corresponding to the resonance lines are identified, through 
their profiles, in the bottom part of Fig. H3 The most intensive line corresponds to bulk mode 
m = 201, having the highest frequency value. The intensity of subsequent resonance lines 
decreases systematically with mode frequency, as can be expected from the growing number 
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FIG. 4: Magnetostatic modes in a rod of length 200a confronted with the spatial distribution of 
the local field H{n) (bold line) depicted along the rod central axis: only absolute values \m + \ mode 
amplitudes are shown here. A striking feature is that the local field line borders the "empty" region 
(where mode amplitudes are suppressed) of bulk- dead modes] the bulk- extended mode frequencies 
lie above the local field range (f2 m > H max ). 



of bulk-extended mode nodes. The resonance intensity nearly vanishes in the comb mode 
frequency region, to revive again and increase in the bulk-dead mode frequency region. Such 
resonance absorption minimum inside a spectrum is very uncommon |8j, and its existence 
may be of practical importance: due to the peak mode density in its vicinity (cf. Fig. 
Ek), such a minimum would indicate a frequency range in which substantial reduction of 
magnetostatic noise could be expected . 

Is there any evidence of the comb mode existence reported in the literature? We believe 
such an evidence to be provided ,, b y Bay, * a, £ in a st ud y based on Reived 
Kerr microscopy, used to detect collective spin wave modes in inhomogeneously magnetized 
NioxiFeo.ig thin-film stripes. The permalloy stripes investigated were 18 nm thick, 2.3 fim 
wide and 1 mm long, and the experiment was performed with an in-plane magnetic field 
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FIG. 5: (a) The discretized magnetostatic mode frequencies (black squares) vs. the mode number, 
m. The bold line represents the spatial distribution of the local field H(n) along the central axis; 
the rod length is 200a. (b) Note the double degeneration of the lower-spectrum {bulk-dead) modes 
with m = 1 — 70 (S/A means, respectively, symmetric-antisymmetric mode). 

oriented along the stripe short axis. A new type of excitation is distinguished by the authors 
among the spin wave modes observed, and named " crossover modes" : though spanning the 
entire stripe, these "new" modes are only detectable at the edges. To explain this feature, the 
authors refer to two different spatial regimes: a small wave- vector region near the edge of the 
stripe, and a large wave- vector region surrounding the stripe center. Although the geometry 
of the samples used in the experiments by Bayer et al. (stripe) is different than that of the 
samples considered in this paper (rod), the properties of the experimentally found "crossover 
modes" strikingly coincide with those of the "comb modes" described here. Moreover, both 
mode types (crossover and comb) have the same origin, appearing as a result of spatial 
inhomogeneity of the local field along the direction of mode formation. The only difference 
lies in the origins of this inhomogeneity: the spatial inhomogeneity of the effective field in 
the samples used by Bayer et al. is due to inhomogeneous magnetization distribution across 
the stripe, while in our study, the effective field inhomogeneity is a consequence of a spatial 
inhomogeneity of the respective demagnetization coefficient (the sample magnetization is 
assumed to be homogeneous). 

Let's get back once more to the resonance spectrum resulting from our calculations and 
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FIG. 6: Resonance spectrum of magnetostatic modes in the rod of the length 200a and the aspect 
ratio 10. Note the existence of the nonabsorbing frequency region corresponding to comb modes. 

depicted in Fig. 6. Note that this spectrum contains a very peculiar segment, corresponding 
to excitation of bulk-dead modes. The peculiarity of this segment lies in the fact that 
the resonance line intensities vary non-monotonically with increasing energy: a rising part 
and a falling one can be distinguished, with a local maximum between them, located in 
the central part of the spectrum. This type of behaviour is very unusual, as resonance 
intensities in multi-peak ferromagnetic resonance spectra tend to be monotonia functions 
of energy. Therefore, we searched the available literature for a similar spectrum segment in 
the results of experimental resonance measurements (performed in rods), and we found an 
old paper by Patton and Schlomann reporting exactly the resonance spectrum segment 
we have been looking for. It appeared as a result of their study quite unusually, and only 
due to the authors' very original idea to distinguish, in their resonance measurements, two 
spectrum components: one associated with the central part of the rod, the other with its 
border parts. 

Patton and Schlomann investigated an axially magnetized cylindrical 0.075 inch (diam- 
eter) 3 inch (length) YIG rod; the central part in which a separate resonance could be 
obtained was 0.4 inch long. The original measurement setup, conceived by the authors, 
allowed a resonance spectrum to be measured first in the whole rod and then, separately, 
only in its central part. The spectrum resulting only from the absorption in the rod borders 
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could be deduced by differentiation from both spectra obtained (Fig. in the quoted study 
Q) A theory of— static modes in axia lly magnetized rods, deveioped by Joseph 
and Schlomann was used as a basis of interpretation of the measurements, but only 
the resonance positions of the magnetostatic modes calculated on this basis were identified. 
The resonance intensities were not interpreted theoretically; the results obtained within the 
present study seem to allow to fill in this gap. We shall focus on the series of resonance lines 
associated with the modes excited in the rod borders as shown in in Fig.4a. Strikingly, 
this resonance pattern is qualitatively identical with the segment of our spectrum (Fig. IHJ) 
identified as associated with bulk-dead modes. Indeed, a deeper sense is hidden behind 
this apparent analogy: since bulk-dead modes contain a dead-amplitude central zone, their 
amplitudes being significant only in the rod borders, obviously, their resonance excitation 
can come from the borders only, exactly as indicated by the Patton and Schlomann's ex- 
periment! Although our theory deals with a square rod, while the Patton and Schlomann's 
experiment was performed on a cylindrical one, we anticipate that this geometrical difference 
is not essential for the above- discussed intensity patterns, since our in-plane sums (jHJ) are 
only slightly modified when a square is replaced by a circle 5]. Thus, anticipating bulk-dead 
modes on a theoretical basis in the present paper, we believe that Patton and Schlomann's 
experiment [10| provides a proof of their existence. 
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